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, Abstract. We show how the renormaUzation group is used to construct a low- 

momentum nucleon-nucleon interaction Mow fc , which unifies all potential models used 
in nuclear structure calculations. T^ow k can be directly applied to the nuclear shell 
model or to nucleonic matter without a G matrix resummation. It is argued that 
Mow k parameterizes a high-order chiral effective field theory two-nucleon force. We 
use cutoff dependence as a tool to assess the error in the truncation of nuclear 
forces to two-nucleon interactions and introduce a low-momentum three-nucleon force, 
which regulates A = 3,A binding energies. The adjusted three-nucleon interaction is 
. perturbative for small cutoffs. In contrast to other precision interactions, the error 

ly-^ I due to missing many-body forces can be estimated, when T^ow k and the corresponding 

■ three-nucleon force are used in nuclear structure calculations and the cutoff is varied. 
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1. Introduction 

There has been much progress over the last five years on improving many-body methods 
3 ■ applicable to nuclei and nucleonic matter. These improvements are most successful 
in different regions of the nuclear chart: the Bloch- Horowitz approach for few-body 



X 



systems, the No-Core Shell Model for light nuclei, the Coupled Cluster Method for the 
intermediate mass region. Density Functional Theory and effective actions for heavy 
nuclei, and the Renormalization Group approach for nucleonic matter. Although these 
microscopic many-body approaches are in principle different methods to diagonalize an 
A-nucleon Hamiltonian, an error is always introduced since it is not possible to include 
up to y4-body forces. Thus, it is important to understand the error of the truncation, 
e.g., to two-nucleon (NN), or two- and three-nucleon (3N) interactions, and to explore 
different choices in the nuclear force starting point. 

When systems are probed at low energies, it is convenient to use low-momentum 
degrees-of-freedom and replace the unresolved short- distance details by something 
simpler, without distorting low-energy observables. As a result, there are an infinite 
number of low-energy potentials corresponding to different resolutions, and one can use 
this freedom constructively to pick a convenient one. In nuclear physics, many issues 
depend on the resolution, e.g., the strength of 3N relative to NN forces, the spin-orbit 
splitting obtained from NN only, or the size of exchange correlations. The change of 
the resolution scale corresponds to changing the cutoff in nuclear forces, and thus this 
freedom is lost if one uses the cutoff as a fit parameter, or cannot vary it substantially. 
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In this Talk, we review results for a "universal" low-momentum NN interaction, 
called V^owfc- This unifies all potential models used in nuclear structure calculations. 
We present Faddeev results for few-body systems and show that the cutoff variation 
of A = 3,4 binding energies is of the same size as results for different precision NN 
interactions, but low- momentum cutoffs are closer to experiment. This demonstrates 
that cutoff independence, and thus model independence, in nuclear physics requires 
consistent 3N forces. After augmenting V^owfc by chiral 3N forces, we find that 3N 
contributions are perturbative for small cutoffs. The set of low-momentum two- and 
three-nucleon interactions can be used in calculations of nuclear structure and reactions 
and we discuss promising directions. Finally, we show that Viowk and G matrix elements 
are quantitatively similar, but Viowk as a potential has a solid theoretical foundation with 
corresponding 3N forces, whereas a G matrix introduces uncontrolled approximations. 

2. Low-momentum nucleon-nucleon interaction 

Conventional precision NN interactions are well-constrained by two-nucleon scattering 
data only for laboratory energies i?iab ^ 350 MeV. As a consequence, details of nuclear 
forces are not constrained for relative momenta k > 2.0 fm~^ or distances r < 0.5 fm. 
However, all these potentials have strong high-momentum components as illustrated by 
the different lines in Fig. ^ This leads to model dependences and technical difficulties in 
many-body applications. Starting from a given potential model Vnn, we have integrated 
out the high-momentum modes above a cutoff A in the sense of the renormalization 
group (RG) PlEj- The resulting low-momentum interaction Viowfc only has momentum 
components below the cutoff and evolves with A so that the low-momentum scattering 
amplitude T{k', k; k"^) (in particular phase shifts and deuteron binding energy) are 
invariant. Thus, in every scattering channel we have 



In order to reproduce the low- momentum T matrix for a given cutoff, V^iowfc is 
renormalized for scattering to intermediate states with p > A. This is achieved be 
resumming high-momentum ladders in an energy-dependent effective interaction, which 
is the solution to the two-body Bloch-Horowitz equation in momentum space with 
projector Q = 6{p — A). The energy dependence can then be recast as momentum 
dependence by using the equations of motion. Both steps are equivalent to the basis 
transformation of Lee-Suzuki. | We note that the largest effect of the renormalization 
is due to the first step of integrating out the high-momentum modes (for laboratory 
energies i^iab ^ 150 MeV the zero-energy Bloch-Horowitz potential describes the phase 

I Note that the RG approach differs from Lee-Suzuki, as we set the Q space biock of the effective 
Hamiitonian (which inciudes aU modei dependences) to zero. 



T{k',k- e) = vMk',k) + ^v ^ 

T{k',k;k') = V,l,{k',k) + -V 




p'^dp, 



p^dp. 



(1) 



(2) 
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Figure 1. Diagonal (left) and off-diagonal (right) momentum-space matrix element 
for V|ow k (symbols) versus relative momentum derived from different high-precision 
potential models for A = 2.1 fm~^. The various bare interactions are given as lines, 
with thick sold or thick dashed lines for the N2L0 (Idaho A) or N3L0 interactions 
respectively. Results are shown for the ^Sg (upper) and ^Si partial wave (lower figures). 



shifts accurately). Both steps are equivalent to integrating an RG equation [T] 
d . w . , _ 2 V,L,{k',A)T\A,k;A') 

For every cutoff l^iowfc defines a new NN potential and a new low-momentum Hamiltonian 

-f^ltwfc = ^ + ^lowfc) (4) 

where the cutoff acts only on the interaction and V^iowfc denotes a (Okubo-) Hermitized 
Viowk (from now on all V^owfc results are for the Hermitian Vio^k and we drop the over- 
line). In many-body applications, H^^^ will lead to different results from T-|-Vnn (since 
unresolved interactions between any high- momentum nucleons are excluded). 

Our main results are shown in Fig. ^ By performing an RG decimation to 
A < 2.1 fm~^, we find that all NN potentials that fit the scattering data and include 
the same long-distance pion physics lead to a "universal" low-momentum interaction 
Viawk mi2|- This holds for all channels and low-momentum off-diagonal matrix elements. 
Note also that in the ^Sq channel, where there is a significant change of Vnn from N2L0 
to N3L0, the Viowfc moves towards the "universal" curve from N2L0 to N3L0. Finally, 
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Figure 2. Evolution of the diagonal Viowfc matrix elements obtained from different 
potentials for cutoffs A 2.0 . . . 4.0 fm^^ versus relative momentum in the ^Si channel. 



we illustrate the collapse in Fig. |21 where we show the evolution of diagonal l^owfc niatrix 
elements from A = 2.0 . . .4.0fm~'^. Further results and details can be found in 

We emphasize that the renormalization of high-momentum modes is theoretically 
and in practice easier in free space, before going to a many-body system.§ Mowfc does not 
require a G matrix resummation, which was introduced because of (mo del- dependent) 
high-momentum modes in nuclear forces. Finally, Viowfc is energy-independent and the 
cutoff is not a parameter (no "magic" value). As we demonstrate in the next Section, 
the cutoff can be used to assess the error of a Hamiltonian truncated to two-body forces. 

3. Cutoff dependence as a tool to assess missing many-body forces 

All NN interactions have a cutoff ( "P-space of QCD" ) and therefore have corresponding 
three- and higher-body forces. Consequently, if one omits the many-body forces, 3N, 
4N,... observables will be cutoff-dependent. Using V\owk all low-energy NN observables 
are cutoff- independent, and therefore, we can assess the effects of the omitted 3N, 4N,... 
forces by varying the cutoff in many-body calculations. In Fig. El we present results 

§ We also note that the RG evolution is very useful for chiral effective field theory (EFT) interactions. 
This is because for lower cutoffs, the phase space for intermediate states is smaller in nuclear structure 
applications. We can start from a chiral EFT interaction with cutoff range, e.g., ~ 500 — TOOMeV to 
include the maximum known long-distance physics, and then run the cutoff down lower. Observables are 
preserved under the RG and higher-order operators are induced automatically, which is more accurate 
and faster than fitting a chiral EFT truncation at the lower cutoff. 




Figure 3. Correlation of the triton and alpha particle binding energies. We contrast 
the Viowk results (left figure) to results for several modern potential models (right 
figure, taken from where plusses denote NN only and diamonds NN with adjusted 
3N forces). Results are given for the T^owfe derived from the Argonne vis or the CD 
Bonn potential. The Tjon line is shown as a linear fit to the NN-only results. The 
thick solid lines in the left figure are interpolations from the Argonne uig/CD Bonn 
results to the respective Mow k for the largest cutoff studied. 



for the and ^He binding energies for a wide range of cutoffs 4]. We find that 
the cutoff variations for all A ^ 1.0 fm^^ are approximately 1 MeV and 5MeV for the 
A = 3 and A = 4 binding energies respectively. The variation should be compared to 
the potential energy, which is (V^owfc) > 30 MeV and > 60 MeV in the two systems. 
The variation gives an estimate of the contribution from many-body forces and shows 
that, while retaining a x^/datum ^ 1, the truncation to a NN potential alone will 
have associated errors of IMeV and 5 MeV for and "^He. We also find that the 
cutoff dependence predicts the linear correlation between binding energies known as the 
Tjon line. Results for reasonable low-momentum cutoffs A ~ 2.0 fm~^ are closer to 
experiment and we observe a slight breaking off the Tjon line. Our results demonstrate 
that 3N forces are inevitable for renormalization and are similar to results obtained in 
the pionless EFT (see proceedings by H.-W. Hammer). Two reference results for the 
Viowk obtained from the Argonne vis potential and A = 1.9 fm~^ are £'(^H) = 8.47 MeV 
and E(^B.e) = 29.19 MeV (for the V^owfc derived from the CD Bonn potential and 
A = 2.1fm-\ we have ^(^H) = 8.49 MeV and E(^He) = 29.20 MeV). 

In the right part of Fig. 01 we also give the potential model dependence in A = 3, 4 
systems. This model dependence is due to the probing of the unconstrained high- 
momentum modes. The high-momentum modes induce three-body correlations, which 
are of short-range. Therefore, at low energies these effects are inseparable from the 
effects due to omitted 3N interactions. Finally, we note that Fujii et al. have reported 
similar results for l^iowfc 1^, but conclude that large cutoffs should be used. However, this 
conclusion is misguided, because it relies on reproducing the few-body binding energies 
of a specific VNN-only model with V^owfc (both without the corresponding 3N forces). 
Since for every cutoff, V^^^^ is a new potential, the few-body binding energies will be 
different, as expected from the right part of Fig. IHl 
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2.5(a) 


33.30 


-40.94 


-2.22 


-0.11 


1.49 
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-90.97 


-11.06 


-0.41 


6.62 


0.12 
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2.5(b) 


33.51 


-41.29 


-2.26 


-1.42 


2.97 


68.03 


-92.86 


-11.22 


-8.67 


16.45 


0.18 


0.74 


3.0(*) 


36.98 


-43.91 


-4.49 


-0.73 


3.67 


78.77 


-99.03 


-22.82 


-2.63 


16.95 


0.23 


0.80 



Table 1. Expectation values of the kinetic energy (T), Viowfc and the different 
3N contributions (long-range 27r-exchange (c-terms), l7r-exchange part (D-term) and 
contact interaction (S-term)) for '^H and ^He. All energies are in MeV and momenta 
are in fm"-'^. (a) and (6) denote two possible solutions for A = 2.5 fm^^ and (*) 
indicates that the ''He fit is approximate, for details see [T. In addition, we give the 
ratio of maximum 3N to Viowfc contribution and an average relative momentum k-^nw 



4. Perturbative low-momentum three-nucleon interaction 

It would be extremely nice to use the RG and calculate a low- momentum NN and 3N 
interaction from different potential models. This however has two problems. First, 
it requires accurate calculations of high-energy 3N scattering wave functions. Second, 
and more severe, with the exception of chiral EFT interactions there is no consistency 
between the NN potential models and their fitted 3N forces. As a consequence, we have 
decided to adjust the leading-order chiral 3N forces to Viowfc for various cutoffs. The 
motivation for this is that at low energies, all phenomenological 3N forces due to meson 
exchanges and high-momentum N, A,... intermediate states collapse to this operator 
form. Therefore, it is reasonable to use the operator form constrained by chiral EFT and 
adjust the coupling constants to Viowfc- Moreover, cutoffs in Viowk and chiral potentials 
are very similar. Both are low-momentum interactions, where only pion exchanges are 
explicitly resolved. When we start from chiral interactions and run the cutoff down 
lower, we find that they also collapse to the "universal" Viowfc, which indicates that 
Mowfc parameterizes a higher-order EFT interaction with sharp-cutoff regularization. 

In chiral EFT, the leading-order 3N interaction enters at N2L0 and consists of 
a long-range 27r-exchange, an intermediate-range l7r-exchange part and a short-range 
contact interaction There are five coupling constants: three low-energy c constants 

in the 2n part, as well as D- and ii^-term couplings in the Ivr and contact term 
respectively. A possible determination of the c constants is through a NN partial wave 
analysis, which includes the long-distance Ivr and 27r physics in the interaction. This has 
been carried out by the Nijmegen group, and we take their values for the c constants [Hj, 
which is most in keeping with our results that Viowfc is strongly constrained by the 
scattering data. Since the c constants parameterize low-energy vrA^ physics, their values 
are independent of the cutoff used to regularize nuclear forces. 

We then adjust the two D- and i?-term couplings to the and ^He binding energies 
for different cutoffs (We note that it may be better to adjust the 3N force to and 
a heavier system, say ^^O, when the latter can be calculated more accurately in the 
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future). For cutoffs A < 2.0 fm^^ we find linear dependences in the fitting, which 
are consistent with a perturbative D- and E-term contribution Ei^B.) = E{Viowk + 
c— terms) + CD{D—teTm) -\-ce{E— term) {cd and Ce are the couphng constants and (. . .) 
denote the matrix elements of the operators). This has been checked explicitly and also 
for the c-terms. For Viowk and all studied cutoffs A < 2.0 fm~\ we find that the adjusted 
3N forces are perturbative by which we mean (\&(^)| VJjnI^I/*^^^) ~ | VsnI^I/'^^^), 
where are exact solutions including up to ra-body forces. We use the operator form 

of the chiral 3N force given in Eq. (2) and Eq. (10) in |7^, multiplied by an exponential 
regulator exp[— ((p^ + 3g^/4)/A^)'^] with the same cutoff value as in Vio^k ij> and q are 
Jacobi momenta). The high power in the exponent yields a behavior similar to a sharp 
cutoff. The two parameters of the 3N force fit to l^iowfc are tabulated for a wide range 
of cutoffs in We note that it is non-trivial that a fit solution of the leading-order 
chiral 3N form with realistic c constants exists. 

Next, we present the different contributions to the triton and alpha particle binding 
energies in Table ^ Assuming that the kinetic energy is due to independent particle 
pairs, we can use this to obtain an average relative momentum fcrms = a/ {k"^) ~ 
mT / (A — 1), when V\owk is used in these systems. Over the range of cutoffs in TablelH we 
find fcrms ~ 0.55 . . . 0.80 fm"^ for ^He (0.50 . . . 0.67 fm"^ for ^H). Although not observable, 
it is reassuring that fcrms ^ A and also intriguing that for all low-momentum cutoffs 
fcrms ~ ^TTi as expected in chiral EFT. We also find that the non-linearities in the fitting 
for larger cutoffs A > 2.5 fm~^ lead to a ratio of the maximum 3N to ^owfc contribution 
K, 0.2. In the chiral counting, 3N contributions are on the order of (Q/A)^ relative to 
the NN force, where Q is a typical momentum in the system. With Q ~ k^^^ ~ m^r, we 
find {Q/KY ~ 0-05 for A ~ 2.0 fm^^ and thus the 0.2 ratio at larger cutoffs is beyond 
this expectation. We take this as an indication that the leading-order chiral 3N force 
is insufficient for larger cutoffs, where more physics is resolved, or that one enters the 
non- linear regime of a limit cycle (see proceedings by A. Nogga). We also observe that 
the c-terms and the ii^-term increase with increasing cutoff and cancel. This is expected 
since the ii^-term renormalizes all divergences in the 3N system. 

It is important to note that the 3N contributions, while perturbative for all cutoffs 
A < 2.0 fm~^, increase by a factor ~ 5 from A = ?> to A = A. This density dependence 
leads to saturation in nuclear matter P|. Finally, for further details on V^owfc in few- 
nucleon systems and the low-momentum 3N force see and proceedings by A. Nogga. 

5. Harmonic-oscillator matrix elements and G matrix comparison 

The advantage of using low-momentum interactions in many-body applications is that 
Viowfc is a soft interaction, without a strong core at short distances. Therefore, ViowA; 
does not couple strongly to high-lying states and can be used in small model spaces. 
This makes it for the first time possible to start directly from a precision NN interaction 
for applications in nuclear structure or reactions. 

The benefit of a lower cutoff in shell model applications is shown in FiglH where 
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Figure 4. Relative harmonic-oscillator matrix elements {n' = ?'|Viow klnl S J) versus 
radial quantum number n for a low-momentum cutoff A — 2.0 fm~^ and a 1.0 GeV 
cutoff A = 5.0 fm~^. Results are shown for the S-wave matrix elements. In both cases, 
Mow*: is obtained from the Argonne vig, potential and Hlu ^ 14MeV. 



we compare S-wave relative harmonic-oscillator matrix elements for a low-momentum 
cutoff A = 2.0 fm"^ and a 1.0 GeV cutoff A = 5.0fm"\ We find that for low- 
momentum interactions the matrix elements decrease quite rapidly and become small 
for \n — n'\ ~ 10. This is not the case for interactions with high-momentum components, 
which require basis states up to ~ 50 shells for convergence. Fig. |3] clearly shows that 
strong high-momentum modes are only poorly represented in a shell model basis. 

In conventional approaches to shell-model effective interactions, the cores are tamed 
by performing a ladder resummation of a Vnn model to obtain a G matrix. However, 
the G matrix resummation introduces an uncontrolled starting-energy dependence and 
requires further approximations in practice. Moreover, there is no theory for the starting 
energy, since the Bloch-Horowitz energy self-consistency is lost when one restricts the 
effective interaction to two-body in an A-hodj system. In Fig.El we compare V^owfc to G 
matrix elements in four major shells. We find that in both T = 0, 1 channels the matrix 
elements are very similar. The biggest differences are on the diagonal T = matrix 
elements. This comes as no surprise since the diagonal matrix elements are related to 
the monopole interaction, responsible for most of the nuclear binding. Here, we expect 
the largest effect of the low-momentum 3N force. A detailed and more quantitative 
comparison will be presented in ^l], where we also study the phenomenology of low- 
momentum interactions and point out where calculations including 3N forces are most 
needed. We note that similar correlations between Viowfc and a G matrix hold for different 
cutoffs and different hw (up to a simple scaling) [TT]. Here, we only want to show that 
Mowfc can be used directly in nuclear structure applications, and emphasize that Viowfc 
has a well-defined theoretical basis with perturbative, low-momentum 3N forces. 



Low-momentum interactions for nuclei 



9 



I • \ • r 
T=0, 4 major shells 
+ diagonal ababJ 




2 - 
- 
-2 - 
-4 - 
-6 - 



-10 -8 -6 -4 -2 2 4 6 8 
G matrix elements, (J)=-80 MeV 



1 ' 1 ' 1 
T=l, 4 major shells 


1 ' 1 ' / 


X diagonal ababJ 




/ , 1 , 1 , 1 


1,1 



-6 -4-2 2 
G matrix elements, co=-80 MeV 



Figure 5. Correlation plots between Mowfc and G matrix elements in 4 major 
shells. The matrix elements {ab\ . . .\cdJT) are in MeV for hui — 14 MeV. We 
also distinguish the diagonal elements related to the monopole interaction and thus 
nuclear binding. Viowk is derived from the Argonne vis potential and the G matrix is 
for Idaho A computed with a rectangular Pauli operator for 4 major shells and starting 
energy —80 MeV ^H] (for other correlation plots see proceedings by B.A. Brown). 



6. Summary and outlook 

In this Talk, we have reviewed advances in constructing low-momentum interactions, 
which can be used directly in many-body apphcations. Our results show that, for low- 
momentum cutoffs, nuclear forces are well constrained and that difficult-to-handle cores 
are not needed to reproduce the NN scattering data. We have shown how the RG can 
be used to construct a model-independent low-momentum interaction, which unifies all 
precision potential models used in nuclear structure calculations. We believe that Viowfc 
is a very useful for nuclear many-body problems for the following reasons: 

(i) It is possible to vary the cutoff in Viowfc (or Viowfc and adjusted 3N interactions) 
over a wide range. This enables one to estimate an error due to omitted many- 
body forces (or omitted higher-order 3N, 4N,... interactions). In this way, it is 
possible to vary the cutoff for interactions with singular pion exchanges and obtain 
an error control in many-body calculations as in the pionless EFT. 

(ii) We have shown that 3N forces are required by renormalization and that adjusted 
low-momentum 3N interactions are perturbative for cutoffs A < 2.0 fm~^. This 
should considerably simplify including 3N forces in nuclear structure applications, 
e.g., shell model interactions, coupled cluster theory or nucleonic matter. 

(iii) Low-momentum interactions bind nuclei in Hartree-Fock ^2], in contrast to all 
other microscopic NN interactions. Consequently, exchange correlations are smaller 
starting from low-momentum forces, and a quantitative derivation of a nuclear 
density functional seems possible [13] ■ Moreover, preliminary results indicate that 
nuclear matter with Viowfc and 3N forces is perturbative for neutron matter the 
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Hartree-Fock equation of state is very reasonable [Tlj . 

potential provides a well-defined starting point for microscopic 
calculations of effective interactions for heavier nuclei in small model spaces. 

Applications of Vio^k that were not discussed in this Talk range from valence particle 
nuclei ^H] to quasiparticle interactions and pairing in neutron matter jTHl El • 

We close with a recommendation of how to start using Vfowfc and some priorities 
for future research. If one wants to use V\owk as a new potential without varying the 
cutoff, we suggest to use cutoffs near A = 2.0 fm~^ (for the I^owA: derived from the 
Argonne f is (CD Bonn) potential the triton binding energy is accidentally reproduced 
for A = 1.9 fm~^ (A = 2.1 fm~^)). We would take these as first cutoff values but 
stress that the 3N force never vanishes. For future investigation, it would be extremely 
promising to study the cutoff variation of nuclear spectra and convergence properties 
using Viovik with 3N forces in the No-Core Shell Model or Coupled Cluster Theory. 
Varying the cutoff will be a powerful tool to provide theoretical error estimates for 
extrapolations towards the drip lines, where one cannot compare to experiment. A 
chart for the spectra of light nuclei with theoretical error bars would be wonderful. 
Finally, more insight on the effects of 3N forces will come from shell model calculations 
using Viovjk with perturbative 3N forces, especially where a two-body G matrix fails. 

It is a pleasure to thank my collaborators Scott Bogner, Gerry Brown, Bengt Friman, 
Dick Furnstahl, Chuck Horowitz, Tom Kuo, Andreas Nogga, Janos Polonyi and Andres 
Zuker for many discussions. This work is supported by the DOE under grant No. DEFG 
0287ER40365 and the NSF under grant No. .nsf-phy 0244822, 
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